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ABSTRACT 

The Einstein Equivalence Principle (EEP) can be probed with astrophysical sources emitting simul¬ 
taneously different types of neutral particles, or particles with varying energies, by testing their time 
of flight through the same gravitational field. Here we use the time delays between correlated photons 
from cosmological transients to constrain the accuracy of the EEP. We take data from two gamma- 
ray bursts as an example, and use, as a lower limit to the theoretical time delays between different 
energies, delays arising from only the gravitational field of our own galaxy. We then show that the 
parameterized post-Newtonian parameter 7 is the same for photons over energy ranges between eV 
and MeV and between MeV and GeV to a part in 10“^, which is at least one order of magnitude better 
than previous limits. Combining this bound on the wavelength dependence of 7 with the absolute 
bound I 7 — 1| < 0.3% from light-deflection measurements at optical (eV) wavelengths, we thus extend 
this absolute bound on 7 to GeV energies. 

Subject headings: 


1. INTRODUCTION 

One statement of the Einstein Equivalence Principle 
(EEP) is that any uncharged test body traveling in 
empty space will follow a trajectory independent of its 
internal structure and composition. For test bodies with 
macroscale masses, the accuracy of the EEP can be mea¬ 
sured in a Newtonian context. However, for test parti¬ 
cles like photons and neutrinos, Newtonian dynamics is 
no longer precise enough. To test the accuracy of the 
EEP at this level, at least the lowest-order deviations 
from Newton’s law of universal gravitation are required, 
which can be conveniently described by the parameter¬ 
ized post-Newtonian (PPN) formalism. For all metric 
theories of gravitation in which all bodies satisfy the 
EEP, the PPN formalism explicitly details the param¬ 
eters in which a general theory of gravity can differ from 
Newtonian gravity. Each theory is specified by a set of 
nume rical coefficients (PPN parameters) see (I Willi [200^ 
I20l4 for a review). 

A gravity theory can be tested through the absolute 
values of PPN parameters. For example, it has been 
pointed out that the time interval required for photons 
to traverse a given distance is longer in the presence of a 
gravitational potential U(r) by 

5t = - [ U{r)dr (1) 

^ J Te 

where denotes the location of the s onrce of parti- 
cle emission and that of observatio n (IShaoirol Il964 
iKrauss fc Tremain3ll988l : iLongollTQSSll . Half of the de¬ 
lay is caused by the warping of space by the gravitational 
field, the other half can be thought of as a result of grav¬ 
itational redshift. Here 7 is one of the PPN parameter, 
representing how much space curvature is produced by 
unit rest mass. Different gravity theories predict differ¬ 


ent values f or 7 . For i nstance, general relativity predicts 
that 7=1 (IWilllll99^ . Many precise methods have been 
devised to test the accuracy of general relativity through 
the measurement of the value of 7 , the most precise be¬ 
ing the measurement of the gravitational deflection of 
light near the Sun and the round-trip travel time delay 
of a n artificial radar signal due to the solar system grav¬ 
ity (I Willi 120061 [Ml : 

• the most precise results from the deflection of light 
is through very-long-baseline radio interferometry 
meas urement, yielding 7 — 1 = (— 0.8 ± 1 . 2 ) x 
10““^ (iLambert fc Le Poncin-Lafittell2009ll2011ll : 

• the most precise results from the travel time de¬ 
lay of a radar signal, which is also the most pre¬ 
cise constraints on the 7 value to date, is from the 
Doppler tracking of the Cas sini spacecraft, yield ing 
7 - 1 = (2.1 ± 2.3) X 10-5 (iBertotti et ahilMl . 

It is worth pointing out that 7 = 1 is not a sufficient 
condition to identify general rela tiyity, since it is not the 
only theory that predicts 7=1 ([Willi 1 1991 . Regardless 
of the absolute value of 7 , all metric theories of gravity 
that incorporate the EEP predict that all test particles 
(photons or neutrinos) must follow identical trajectories 
and undergo the same time delay, even though the delay 
time may vary between different theories. In other words, 
as long as the EEP is valid, all metric theories predict 
7 = 7 i = 72 (and 7 need not even be 1 ), where 1 and 
2 represent two different test particles. In this case, the 
accuracy of the EEP can be characterized by limits on 
the differences of 7 value for different species of particles, 
or the same specie of particle with varying energies. 

It has been proposed that the close coincidence in time 
of arrival of the photon and neutrino bursts from the su- 
pernova 1987A provides a strong test of the EEP ([Longol 
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Il988t iKrauss &: Tremaindll988l) . Compari ng: the arriva l 
time delay of the photons and neutrinos, iLongol (Il988f ) 
proposed that the 7 values for photon (eV) and neutrino 
(MeV) are identical to an accuracy of approximately 
0.37%. On the other hand, the very small dispersion 
in arrival times observed for neutrinos of different ener¬ 
gies (7.5 — 40 MeV) can be used to set a limit on the 
energy dependence of "fu to the level of 1.6 x 10 “®. 

Supernova 1987A is an extragalactic source, but non- 
cosmological. In this work, we propose that the time de¬ 
lays between correlated photons from cosmological tran¬ 
sients, such as Gamma-Ray Bursts (GRBs), can also be 
used to constrain the EEP with similar approaches. In 
the following, we first discuss the general picture of using 
cosmic transients to constrain the EEP. Taking GRBs as 
an example, we then show the ability of this method for 
improving the constraints, and discuss the uncertainties 
involved. 

2. GENERAL DESCRIPTION FOR DERIVING EEP 
CONSTRAINTS WITH COSMIC TRANSIENTS 

Gosmic transient sources often show some time delays 
between different energy bands. The main terms that 
might contribute to this time delay may be expressed as 
follows: 


Atobs — Afint + Atpiv + Atgpe + Atg^a (2) 

In the following, we discuss each component in turn. 

• Afint is the intrinsic (astrophysical) time delay be¬ 
tween two test photons, and this has the largest 
uncertainty since it is impossible to determine it 
based on the observational data alone. Appropri¬ 
ate assumptions for the radiation mechanism are 
usually required to assign a chronological order for 
the different energy photons, but it is still hard to 
estimate the exact value of Afjnt. 


• AtLiv represents the potential time delay due to 
Lorentz invariance violation (LIV) via an energy- 
dependent velocity of light. Some descriptions of 
it are based on quantum fluctuations in the back¬ 
ground space-time metric causing a non-trivial dis¬ 
persion in vacuo for lig:ht propagating through thi s 
space-time foam (|Amelino-Gamelia et al.l Il998f) . 
Taking into account the effect on the propagation of 
photons due to the expansion of the Universe, even 
a tiny variation in ph oton speed may lead to an ob¬ 
serva ble time delay (lEllis et al.l 120061 : lAbdo et'aP 

ISMi, 


(1 -b n) E'l - r (1 -b zTdz' 

~ 2i7o Jo Hz') 


( 3 ) 


which is the n-th (first non-vanishing term) of a 
series expansion in an effective field theory formu¬ 
lation of quantum gravity including Lorentz invari¬ 
ance violating terms. Here Ei and E 2 represent 
the energies for two test photons {Ei > E2), Eqq 
is the effective energy scale of quantum gravity, 
h{z) = H{z)/Ho, with Hq being the present Hub¬ 
ble expansion rate and H{z) the Hubble parame¬ 
ter at different redshifts. For a limiting Eqq equal 
to the Planck energy, the recent evidence disfavors 
the existence of linear (order n = 1) LIV terms, i.e. 


energy-dependent departures of t he speed of pho¬ 
tons from the usu a l speed of light (fEllis et aLll2006l 
lAbdo etld] l2009t iPan et aH 1201^ Substituting 
the Plank energy for Eqq it is easy to show that 
the higher order terms of the expansion (n > 2 ), 
even if they exist, are negligible for the purposes 
of this work We thus ignore this term in the 
following analysis. 


• Atspe represents the potential time delay due to 
special-relativistic effects in the case where the pho¬ 
tons have a rest mass which is non-zero, which 
reads 


/ I 1 (1 + z')dz' 

2Ho [Ef El) h{z') 


( 4 ) 


where rUph is the photon rest mass. Modern exper¬ 
iments have provided the upper limit s for the pho¬ 
ton r est mass as ruph < 10“^® eV/c^ (|Amsler et al.l 
I 2 OO 8 D . In this case, AGpe is negligible even when 
the energy of test photons are lower than radio 
band. 


• Atgra corresponds to the difference in arrival time 
of two photons of energy Ei and E 2 , caused by 
the gravitational potential U{r) integrated from 
the emission site to the Earth, which can be es¬ 
timated with equation [TJ as long as the expression 
of U{r) is known, 

Atgra = 3 -^ f U{r)dr (5) 

^ Jto 

where c = 3 x 10^° cm s“^ is the energy- 

independent speed for massless particles. In gen¬ 
eral, U(r) can be divided into three parts: the grav¬ 
itational potential of our galaxy t/mw(?') would be 
dominant out to a certain distance, say to the edge 
of the local group (^ 1 Mpc). Based on the cosmo¬ 
logical principle and current CMB results, the grav¬ 
itational potential should gradually approach a flat 
background intergalactic potential Uicir) between 
our galaxy and the host galaxy (or cluster) of the 
transient. This assumption may be valid at least 
for sources with z < 2, for which th e lensing prob¬ 
abilit y is significant only to < 1% (|Wiklind et al.l 
l2013f) . The final contribution is from the gravity 
well of the host galaxy (or the relevant galaxy clus¬ 
ter) t/host(7’)- 

Leaving out the negligible components, equation!^ is 

Atobs = Atint + Atgra ( 6 ) 

For a specific case where high energy photons arrive later 
than the low energy photons (Afobs > 0 ), if we know that 
Afint > 0 , we will have 

Atobs > ^3 J U(r)dr (7) 

The potential function for UiG{r) and t/host(?’) is hard to 
model, but it is very likely that the effect of these two 

^ In general, one can find that as long as Eqq^^iIE i ^ 
the n-th expansion of AtLjv becomes negligible for the purposes 
of this work (e.g., much smaller than the order of a second). 
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terms is much larger than if we simply assumed that the 
potential is just ?7mw(r) extended to the distance of the 
transient, i.e., 

Atobs > 3^^ f U^-w{r)dr (8) 

^ «/To 

Although the gravitational potential of the Milky 
Way at large distances is still not well known, 
iKrauss fc Tremaini (|1988ll examined two popular poten¬ 
tial models, e.g., the Keplerian potential U{r) = —GM/r 
and the isothermal potential U{r) = [ln(r/rniax) — 1], 
for r < r^nax, where Vc is the circular speed at the 
solar radius and the potential becomes Keplerian for 
r > Cmax = 100 kpc. They proposed that different mod¬ 
els for Univj{r) do not have a strong influence on the re¬ 
sults for testing EEP. In this work, we adopt the Keple¬ 
rian potential for our galaxy. Thus, for the purposes of 
obtaining a lower limit, we can simply extend this galac¬ 
tic potential out to cosmic scales to bracket from below 
the potential function of Uioir) and Uhost{r)- We thus 
have 

7 .- 7 .<AU.(^)"ln-(^) (9) 

where G = 6.68 x 10“® erg cm g~^ is the gravitational 
constant, ~ 6 x is the mass of our galaxy 

(|McMillanll2011l : iKafle et al.ll2012f) , d is the distance from 
the transient to Earth, and b is the impact parameter of 
the light rays relative to the center of the galaxy. For a 
cosmic source in the direction (R.A. = /3s, Dec. = 6s), 
the impact parameter can be estimated as 

b = rG\/l — (sin^s sin^c -|- cos 6s cos 6 g cos(/3s — /3g))^ 

( 10 ) 

where tq = 8.3 kpc is the distance from the Sun to 
the galaxy center, and (/3 g = 17*'45™40.04®, 6o = 
—29°00'28.1") are the coordinates of the galaxy center in 
the e quatorial coordinate system fJ2000: [Gillessen et al.l 
1200911 . Note that a more accurate formulation for Umwir) 
may provide higher order corrections to equation!^ but 
the correction would be absorbed into the approxima¬ 
tion obtained from replacing UiQ{r) and C/host(?') with 
the extension of 

3. EEP CONSTRAINTS WITH GRBS 

Gamma-ray bursts are the most extreme explosive 
events in the universe. They were initially observed as 
short, intense, and non-repeating flashes of ^ MeV 7 - 
rays. Later on, it was found that the initial burst is 
usually followed by a longer-lived broadband emission 
(X-ray, ultraviolet, optical, infrared, microwave and ra¬ 
dio), showing that GRBs are actually multi-wavelength 
transients, instead of simply bursts of 7 -rays. The initial 
MeV 7 -ray emission is often called the “prompt emission, 
while the longer wavelength emission is referred to as the 
“afterglow emission. It is worth noticing that in some 
bursts, the longer wavelength emission, e.g. optical emis¬ 
sion, sometimes already starts during the prompt emis¬ 
sion phase. The duration of the prompt emission is rela¬ 
tively short, ranging from ^ 0.1 s to ~ 1000 s. Morever, 
the histogram of the duration shows a bimodal structure 
separated at 2 s, where bursts with duration less than 
2 s are classified as short-GRBs and those that last for 


more than 2 s are called long- GRBs. (|Kumar fc Zhan j 
I 2 OI 4 for a review). 

GRBs are interesting cosmic transients for obtaining 
EEP constraints for the following reasons: 

• Even the 7 -ray photons are usually observed in dif¬ 
ferent energy bands, for instance, the two instru¬ 
ments onbo ard Fermi, the Gamm a-ray Burst Mon¬ 
itor (GBM; iMeeean et al.l (l2009fll and the Large 
Area Telescope fLAT: I Atwood et al.l (|2009f B. pro¬ 
vide unprecedented spectral coverage for seven or¬ 
ders of magnitude in energy (from ^8 keV to ^300 
GeV). 

• The observed 7 -ray lightcurves in each energy 
band often have violent variability, wit h minimum 
times cales as small as milliseconds (iBhat et al.l 
mH). It is easy to use these sharp features to 
identify the arrival time lag between different en¬ 
ergy photons. 

• Thanks to the prompt slewing capabi li ty of the X- 
Ray Telescope (XRT; iBurrows et al.l (12005111 and 
the UV-Optical Telescope (UVOT; Roming et al. 

S ) onboard the Swift satellite ( Gehrels et al. 

longer wavelength photons have been ob¬ 
served at very early times after the GRB trigger, 
even during the prompt emission phase. It is thus 
possible to determine the arrival lag between 7 -ray 
photons and low energy photons, such as optical 
photons. 

After four decades of observation, thousands of GRBs 
have been studied in significant detail, making it pos¬ 
sible to find some good examples that can be used to 
constrain the EEP. In this work, we take two examples, 
GRB 090510 and GRB 0803I9B. 


3.1. GRB 090510 

GRB 090510 Hrst triggered the Fermi-GBM on a pre¬ 
cursor at To = 00 : 22 : 59.97 UT, May 10th, 2009, 
and the main emission episode in the 8 keV — 40 MeV 
energy range started ^ 0.5 s after Tq, lasting up to 
^ 1 s. Then, 0.43 s after the GBM trigger, it was 
detected and precisely localized by the Swift satellite, 
with coordina tes (J2000) R.A. = 22''14“12.47®, Dec. = 
—26°35'00.4" (|Hoversten et al.l[200911 . The emission ob¬ 
served by the Fermi- LAT starte d 0.65 s after the 
GBM trigger an d lasted ~ 200 s (|Guiriec et al.l 120091 : 
lOhno fc Pelas 'iJ [200^ . The Tgo of the burst is 0.3 s, 
placing it in the short GRB category. Subsequent ob¬ 
servations dAterminedJts_^edslnfta^ = 0.903 ± 0.003 
(iRau et al.ll2009l: iMcBreen et al.l[2010ll . 

lAbdo et al.l ( 2009f ) analyzed in detail its prompt light 
curve over the different energy bands observed by Fermi. 
They found that the total GBM light curve (8 keV — 
40 MeV) consists of seven main pulses, where the Hrst 
one is a dim short spike. A single 31-GeV photon was 
detected by LAT at 0.83 s after the GBM trigger, which 
coincides in time with the last of the seven GBM pulses. 
The directional and temporal coincidence of this photon 
with the MeV radiation of GRB 090510 is very si gnifi- 
cant, at the > 5a confidence level (|Abdo et al.ll2009r) . No 
high-energy (GeV) photon has ever been detected before 
the onset of the low-energy (MeV) emission in a GRB, 
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and this 31-GeV photon similarly appears likely to have 
been emitted after the onset of GRB 090510. Taking as 
our nominal time delay At ~ 0.83 s, we thus have from 
eq.® for this burst 

7GeV — 7MeV < 2 X 10“® (11) 

3.2. GRB 080319B 

GRB 080319B was detected by Swift at Tq = 
06 : 12 : 49 UT, March 19th, 2008, with coordi- 
nates (J2000) R.A. = 14>^31”^40.7", Dec. = +36°18'14.7" 
(jRacusin et al.l[^08a^ . It is a long burst with Tgn > 50s 
and located dX z = 0.937 (|Vreeswiik et al.l [200^ . At 
the trigger time of this burst, the wide-held robotic op¬ 
tical telescope “Pi of the Sky” (|Cwiok et al.ll2007tl and 
the wide- held robotic instrument Telescopic Ottimiz- 
zato pe r la Ricerca dei Tran sient! Ottici Rapid! (TOR- 
TOR A; [MohnarreF^ (|2006[ 1'1 both serendipitously had 
the GRB within their helds of view (as they were both 
already observing GRB 080319A). A bright optical tran¬ 
sient was observed starting at 2.75 ± 5s and peaking 
around 18s after the Swift trigger. The peak V-band 
magnitude reache d 5.3, corresponding to a hux den - 
sity of ~28.7 Jy (iKarpov et al.l 120081 : IZou et al.ll2009ll . 
TORTORA measured the brightest portion of the opti¬ 
cal hash with high time resolution, distinguishing three 
separate peaks which enabled doing detailed compar- 
isons between the p romp t optical and 7 -ray emissions. 
iKumar fc NaravanI (j2009t l calculated the optical- 7 -ray 
correlation function for the entire prompt emission, and 
they found that for the contemporaneous rise at 0 — 18 
s and fall at 43 — 60 s, the correlation coefficient has a 
maximum at lags between 0 and 5 s , while for the middle 
part of the prompt emission, the correlation function still 
has a peak at lags between — 1 to 3 s, but is less promi¬ 
nent. To be conservative, we adopt the largest value 5 s 
as the arrival lag At between the MeV photons and the 
optical photons. From eq. we thus have for this burst 

7eV — 7MeV < 1-2 X 10“^ (12) 

4. DISCUSSION 

In the last two sections, we have proposed a practical 
method for using cosmic transients for testing the energy 
dependence of the 7 parameter quantifying departures 
from the EEP, and we have taken GRBs as an example 
to illustrate the capabilities of this method. In the fol¬ 
lowing, we discuss some of the uncertainties involved in 
this method and the caveats on our constraint results. 

The largest uncertainty is from the intrinsic time de¬ 
lay between two test photons, Afint. Since we are only 
seeking upper limits, in the calculations we neglected the 
Atint term based on the assumption that Afint has the 
same sign as Atobs- More specifically, here we assumed 
that the GeV (optical) photons are intrinsically emitted 
later than (or at most simultaneously with) the MeV pho¬ 
tons in GRB 090510 (GRB 080319B). Erom the observa¬ 
tional point of view, such an assumption is empirically 
justified, since we have never observed any GRBs with 
GeV or optical photons arriving earlier than the MeV 
photons. In standard GRB models, the MeV photons are 
assumed to arise in internal shock regions via synchrotron 
emission, while GeV photons are usually believed to orig¬ 
inate from the inverse Compton up-scattering of MeV 


photons, and the optical photons are usually assumed 
to arise in the external shock region (which is outer 
and occurs later than the internal shock) also via syn¬ 
chrotron emission. In other words, from the theoretical 
point of view, our assumption about the intrinsic order of 
GeV or optical photons with MeV photons is also stan¬ 
dard practice. However, if in reality the GeV or optical 
photons are emitted earlier than the MeV photons, and 
|Atint| |Atobs|, our constraint results would be invali¬ 
dated. 

The second uncertainty comes from the determination 
of Atobs- For GRB 090510, we chose Atobs as the time 
delay between the 3I-GeV photon and the trigger pho¬ 
tons for the Swift satellite. Our results could be more 
stringent if in reality the 31-GeV photon were correlated 
with the other GBM (MeV) pulses, for example the sev- 
enth pulse which c oincides in time with the GeV photon 
(|Abdo et al.l 1200911 . However, in principle, it would be 
conceivable that the GeV photon is correlated with some 
other MeV photons emitted earlier but which failed to 
trigger the Swift satellite. This would be similar to a case 
where intrinsically this GeV photon is emitted earlier 
than the MeV photons, which could potentially invali¬ 
date our results. For GRB 0803I9B, we adopted the time 
delay results from the optical- 7 -ray correlation function. 
These results should be reliable, since the TORTORA 
telescope provided high time resolution measurements 
for the brightest portion of the optical flash. Further¬ 
more, it is worth noting that our adopted delay time is 
even larger than the time delay between the onset times 
of the optical observations and the MeV trigger. 

The total mass and the exact gravitational potential 
function of our galaxy is not well known. More accurate 
values for Mmw or for the function C/mw(r) could improve 
our constraints, but the correction should be limited to 
less than an order of magnitude. 

In the calculations presented we neglect the time delay 
terms caused by the host galaxy and the intergalactic 
background gravitational potential. In principle, these 
terms could be much larger than Atgra.m-w- With fu¬ 
ture observations, and with a better understanding of 
the potential functions for these terms, our result could 
be improved by orders of magnitude. 

5. CONCLUSION 

In principle, the accuracy of the EEP may be charac¬ 
terized by limits on the differences between the 7 val¬ 
ues (one of the PPN parameters) for different species of 
particles, or for the same specie of particle at different 
energies. In this work, we propose that the time de¬ 
lays between correlated photons from cosmological tran¬ 
sients, such as GRBs, can be used to constrain the ac¬ 
curacy of the EEP. Using the data from two observed 
GRBs as an example, and assuming that the observed 
time delays between photons of different energies are 
caused dominantly by the gravitational effects of our 
galaxy, an assumption which we justified in § 2 , for the 
GeV and MeV photons of one burst we give an up¬ 
per limit of 7GeV “ 7MeV < 2 X 10“®; and for the 
optical and MeV photons form another burst we find 
7eV — 7MeV < 1-2 X 10“"^. 

Previous to this present work, the most accurate con¬ 
straints on 7 differences were based on the photons and 
neutrinos from Supernova 1987A, giving an upper limit 
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of 0.37% for optical photons and MeV neutrinos, and an 
upper limit of 1.6 x 10 “® for two neutrinos, whose energy 
difference was of order ^ MeV. If our present interpre¬ 
tation is correct, our results increase the accuracy level 
of such EEP constraints by at least one order of magni¬ 
tude, to ^ 10 “^, while extending the energy range tested 
to the GeV-MeV and MeV-eV range. 

A new type of transients, calle d Fast Radio 
Burs ts, have lately attracted attention (iThornton et al.l 
l2013ll. Although the ir physical origin is still debated 
( Kulkarni et al.ll2014ll . if in the future they are proven 
to be cosmic transients, their simple sharp features make 
them attractive candidates for constraining the EEP, 
which could further extend the tested energy range to 
the radio band with high accuracy. 

It is worth pointing out that high energy neutrinos 
have long been proposed to be associated with GRBs, 
although a dedicated search of high-energy neutrinos co¬ 
incident with electrom agnetically detected GRBs has so 
far le d to null results (|Abbasi et al.ll 201 ^ lAartsen et al.l 
l2015ll . If in the future GRBs with both high energy pho¬ 
tons and neutrinos are detected, better constraints on 
EEP might be achieved. On the other hand, neutrino 
detectors such as IceGube have reached the sensitivity 
to detect high-energy neutrinos of astrophysical origin 
for the first time. The three-year discovery catalog of 
IceGube contains a total of 37 neutrino candidate events 
with deposited ener gies ranging from 30 to 2000 TeV 
(|Aartsen et al.l 12014) . Although the origin of these neu¬ 
trinos is still debated, some of them seem to have both 
spatial and tempor al correlations , such as events 24 and 
25 in the catalog (lAartsen et al.l 120141) . In the future, 
if the origin of these neutrinos is better understood, or 
the association between different detected events is con¬ 
firmed, it will be essential to use these high energy cosmic 
neutrinos for obtaining better constraints on the accu¬ 


racy of the EEP. 

The most precise previous constraints on the absolute 
7 value of photons isy — 1 ^ 10 “®, which comes from 
radio photons. Our relative limits are three orders of 
magnitude lower than this, focusing on the value of the 
difference between various energies, rather than the ab¬ 
solute value. On the other hand, the Shapiro time-delay 
measurements using the Hipparcos optical astrometry 
satellite yielded an agreement with Ge neral Relativity of 
7 eV — 1 ^ 0.3% (|Froeschle et aLlHOQTl) . Based on our re¬ 
sults, we can predict that the constraint on the absolute 
7 values MeV or GeV photons, which from the sources 
we investigated here differ from that of their correspond¬ 
ing optical photons by < 10”’’, should agree with GR at 
least to this same < 0.3% level of accuracy. Thus, we 
have extended the validity of this absolute value bound 
on 7 at this level from the optical to the GeV range, 
and have constrained the relative variation of 7 between 
optical and GeV to the <10”'’ range. 
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